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Spin squeezing and precision probing with light and samples of atoms in the gaussian 

approximation 
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We consider an ensemble of trapped atoms interacting with a continuous wave laser field. For 
sufficiently polarized atoms and for a polarized light field, we may approximate the non-classical 
components of the collective spin angular momentum operator for the atoms and the Stokes vectors 
of the field by effective position and momentum variables for which we assume a gaussian state. 
Within this approximation, we present a theory for the squeezing of the atomic spin by polarization 
rotation measurements on the probe light. We derive analytical expressions for the squeezing with 
and without inclusion of the noise effects introduced by atomic decay and by photon absorption. 
The theory is readily adapted to the case of inhomogeneous light-atom coupling [A. Kuzmich and 
T.A.B. Kennedy, Phys. Rev. Lett. 92, 030407 (2004)]. As a special case, we show how to formulate 
the theory for an optically thick sample by slicing the gas into pieces each having only small photon 
absorption probability. Our analysis of a realistic probing and measurement scheme shows that 
it is the maximally squeezed component of the atomic gas that determines the accuracy of the 
measurement. 
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I. INTRODUCTION 

With spin squeezed atomic ensembles, i.e., samples 
where the variance of one of the angular momentum 
(spin) components is reduced compared with the coher- 
ent state value, one has the possibility to measure cer- 
tain atomic and/or classical parameters beyond the pre- 
cision set by the standard quantum noise. Recent ex- 
amples where this possibility was exploited include stud- 
ies on magnetometry with collective atomic spins [l|,Q|. 
The central feature in those works is the entanglement 
of collective continuous light-atom variables. This en- 
tanglement can be created by the free-space interaction 
between a trapped polarized atomic sample and an ap- 
propriately polarized propagating laser beam with pho- 
ton energy adjusted to the energy spacing between the 
atomic energy levels [j,|j|. The probing of the atomic en- 
semble with the light field squeezes the atomic observable 
(the atomic spin) and enables an improved measurement, 
e.g., of a magnetic field. The underlying squeezing of the 
collective atomic spin variable was dealt with in a series 
of papers (see, e.g., Refs. ^5, 6, 7, 8, 9, 10, IJ, and 
references therein) including investigations of quantum 
non-demolition feedback schemes |1(1 112| , and a study of 
the case of inhomogeneous light-atom coupling [llll. In 
the present work, we follow the lines of Refs. |all.^ll4|. 
and investigate the spin-squeezing of continuous variable 
quantum systems in the approximation where the atomic 
and photonic degrees of freedom are described by a gaus- 
sian state. To this end we will use that the gaussian state 
is fully characterized by its expectation value vector and 
its covariance matrix and we will use that explicit formu- 
lae exist for the time evolution of the system and for the 
quantum state reduction under measurements, see, e.g., 
Refs. \WL iwL Il7l | and references therein. In particular, 



the fact that the measurements are explicitly accounted 
for in the gaussian approximation is a strength of the 
present theory. 

In the development of the theory, we shall consider a 
continuous wave (cw) beam of light passing through a 
cloud of trapped atoms. In the Schrodinger picture we 
have an explicit update formula for the quantum state 
conditioned on the outcome of measurements carried out 
on a quantum system, but a light beam is a multi-mode 
field with an infinite dimensional Hilbert space, in which 
a complete description of the quantum state is normally 
prohibitively complicated. The quantum mechanical de- 
scription of cw optical fields is often formulated in terms 
of temporal correlation functions or the noise power spec- 
trum of field operators in the Heisenberg picture, which 
is, however, not a convenient formulation, when the field 
is being monitored continuously in time. When we re- 
strict ourselves to gaussian states, however, it is possible 
to describe the field in the Schrodinger picture and to 
dynamically evolve the combined quantum state of the 
interacting light field and atomic system. 

The paper is organized as follows. In Sec.^ we derive 
the Hamiltonian for the collective atom-light coupling. 
In Sec. mil we describe dynamics and measurements in 
the gaussian approximation and provide update formulae 
for the covariance matrix and for the expectation value 
vector. In Sec. IIVL we present fully analytical results for 
spin-squeezing of an atomic gas for a homogeneous light- 
atom coupling, and small photon absorption probability 
and atomic decay rate. In Sec. we describe how to 
handle the case of inhomogeneous light-atom coupling. 
In Sec. IV Al we treat the case of an optically thin gas, 
i.e., small photon absorption, and we obtain analytical 
results. In Sec. IV BL we investigate the case of an opti- 
cally thick gas. In Sec. IVII we show that the maximally 



squeezed component of the gas will set the limit for the 
precision in a given measurement. Sec. I VIII brief! v sum- 
marizes the results and concludes the paper. 



II. COLLECTIVE ATOM-LIGHT COUPLING 

To describe the atom-Hght coupHng, we imagine that 
the beam is split up into short segments of duration r and 
corresponding length L = ct. These beam segments are 
chosen so short, that the field in a single segment can be 
treated as a single mode, and that the state of an atom 
interacting with the field does not change appreciably 
during time r, so that the evolution of the atomic system 
is obtained by sequential interaction with light segments. 
Since we are interested in modeUing a cw coherent beam 
with constant intensity, we assume a mode function for 
each segment of the field which is constant on a length L 
and within the transverse area A, i.e., the quantization 
of the field energy LAeqE^ = Nphhu yields the relation 
between the electric field amplitude and the photon num- 
ber in the segment of the field, E — \/N^Jj 



hu 



LAeo 



. In the 



scheme for spin squeezing, we consider a fight beam lin- 
early polarized along the x direction and propagating in 
the y direction. The polarization can be decomposed in 
two polarization components with opposite circular po- 
larization with respect to the quantization axis z. These 
two components interact differently with the atoms be- 
cause of the selection rules of the optical dipole transi- 
tion. Imagine atoms with a ground {\g)) and an excited 
(|e)) state with J = 1/2, interacting with the cr+ and 
(T"~ components of the light field on the I5-1/2) ^^ lei/2) 
and |(?i/2) ^^ |e_i/2) transitions, respectively. The inter- 
action Hamiltonian between a collection of iVat atoms, 
enumerated with the index i and the two quantized fields 
thus writes 

H = ^{hga+\ei/2,i){g-i/2,i\ + h.c. 

-f /iga- |e_i/2^i) (51/2,,; I + ^-C-) : (1) 

with fig — ~dEo, d the atomic dipole m oment on the 
relevant transition, and Eq = ^Jhw/LAsQ the 'field per 
photon', identified above. We assume that the fields are 
frequency detuned by an amount A with respect to the 
atomic resonance. In the Hmit where g^/N^ <C A the 
atoms are not excited by the fields, and the dynamics is 
entirely associated with the light-induced energy shifts 
of the ground states. Adiabatic elimination of the upper 
states then leads to the effective Hamiltonian 



H = "^-^ [ala+\g^l/2,^){g^l/2,^\ + 

i=l 

+ a-a-l5i/2,,:)(5i/2,j|j , 



(2) 



which applies for the duration r for which the field over- 
laps the atomic system. The photon field is suitably 



described by a Stokes vector formalism, with a macro- 
scopic value of the component {Sx) = hNph/2, and 
where the Sz operator yields the difference between the 
number of photons with the two circular polarizations, 
Sz = h{a\_a+ — a_a_)/2, and Sy yields the difference be- 
tween the number of photons polarized at 45°and 135°, 
with respect to the z axis, respectively. The Stokes vector 
components obey the commutator relations of a fictitious 
spin, and the associated quantum mechanical uncertainty 
relation on Sy and Sz, Var(5'j,)Var(S'^) = \{hSx)\'^/4:, are 
in precise correspondence with the binomial distribution 
of the Hnearly polarized photons onto the other sets of 
orthogonal polarization directions. We introduce the ef- 
fective cartesian coordinates 



(Xph^Pph) 



Oy 



Sz 



VWsJl^ VWsJl 



(3) 



with the standard commutator [a;ph,Pph] — i and result- 
ing uncertainty relation, which is minimized in the initial 
state, implying that this state is a gaussian state, i.e., its 
Wigner function is a gaussian function of the phase space 
coordinates. 

The atomic ensemble is initially prepared with all N^^^ 
atoms in a superposition {\g-i/2) + |5i/2))/'s/2 of the 
two ground states with respect to the quantization axis 
z, i.e., the total state of the atoms is initially given by 

((I5-1/2) + l5i/2))/\/2) "*• In this state, the system of 
two-level atoms is described by a collective spin vector, 
where the component along the x-direction attains the 
macroscopic value {Jx) = hNatf^, and where the collec- 
tive spin along the z-axis, Jz, represents the population 
difference of the |5±i/2) states. As for the photons, the 
quantum mechanical uncertainty relation for the collec- 
tive spin components of the atomic state corresponds ex- 
actly to the binomial distribution of the atoms on the two 
ground states, and also here it is convenient to introduce 
cartesian coordinates 



(a; at, Pat) = 



J. 



( Jy 



(4) 



for which the initial state is a minimum uncertainty gaus- 
sian state. 

The Hamiltonian ^ can be rewritten in terms of the 
effective atomic and field variables. First, we note that 

J2f=i l5Tl/2,^)(ffTl/2,^l = A^at/2 ± Jj^i and that ala± = 
<i>r/2 ± Sz/fi, where $ is the photon fiux. We then insert 
these expressions in Eq. ^, leave out a constant energy 
shift and obtain the effective interaction Hamiltonian 



Ht = hUrPi^tPph- 



(5) 



We display the product of H and r, to expose the effect of 
the interaction with the whole segment, and we introduce 
the effective coupling 'constant' 



Kr =2- 



\{J.)\\{S.)\ 



(6) 



The free-space coupling constant of light and atoms is 
small, and the coarse grained description will be per- 
fectly valid even for the macroscopic values of A^ph = $t 
required by our treatment. The Hamiltonian in Eq. Q 
correlates the atoms and the light fields. It is bilinear in 
the canonical variables, and hence preserves the gaussian 
character of the joint state of the system [T^. We have 
emphasized the convenience of using gaussian states, be- 
cause their Schrodinger picture representation is very effi- 
cient and compact. Now, given that every segment of the 
optical beam becomes correlated with the atomic sample, 
as a function of time, the joint state of the atom and field 
has to be specified by a larger and larger number of mean 
values and second order moments. If no further interac- 
tions take place between quantum systems and the light 
after the interaction with the atoms, there is no need to 
keep track of the state of the total system. In practice, 
either the transmitted light may simply disappear or it 
may be registered in a detection process. In the former 
case, the relevant description of the remaining system is 
obtained by a partial trace over the field state, which 
produces a new gaussian state of the atoms. We are in- 
terested in the case, where the polarization rotation of 
the field is registered, i.e., the observable Xph is mea- 
sured. The effect of measuring one of the components 
in a multi-variable gaussian state is effectively to pro- 
duce a new gaussian state of the remaining variables as 
discussed in detail in Sec. IIIII 



III. 



DYNAMICS AND MEASUREMENTS IN 
THE GAUSSIAN APPROXIMATION 
INCLUDING NOISE 



Having established the fact that the quantum state of 
the atoms is at all times described as a gaussian state, we 
shall set up the precise formalism. For the column vector 
of the four variables y = {x^t, Pat, Xph, Pph)'^ describing 
the atoms and a single segment of the light beam, the 
Heisenberg equations of motion yield 



y{t + t) = Sry{t) 
with the transformation matrix 



S. = 



1 K^ 

10 

Kr 1 

1 



(7) 



(8) 



From Eq. Q and the definition of the covariance matrix 
7ij = 2Re ((y, - {yi))iyj - (yj))) [llli3, we directly ver- 
ify that 7 transforms as 



jit + T) = Srjit)Sr' 



(9) 



due to the atom-light interaction. 

In the probing process there is a small probability that 
the excited state levels which were adiabatically elimi- 
nated from the interaction Hamiltonian of Eq. ^ will 



be populated. If this happens, the subsequent decay to 
one of the two m^ = ±1/2 ground states occurs with the 
rate 






(10) 



where F is the atomic decay width and a = A^/(27r) is 
the resonant photon absorption cross-section. The con- 
sequence of the decay is a loss of spin polarization since 
a detection of the fiuorescence photons in principle could 
tell to which ground state the atom decayed. If every 
atom has a probability rj-r = r]T to decay in time r with 
equal probability into the two ground states, the collec- 
tive mean spin vector is reduced by the corresponding fac- 
tor (J) -^ («^)(1 — Vt)- When the classical x-component 
is reduced this leads to a reduction with time of the cou- 
pling strength Kr 1-^ Kt-V^ ~ Vt which was also discussed 
in Refs. 0, QjllJI- Simultaneously, every photon on its 
way through the atomic gas has a probability for being 
absorbed fl4 | 



N,,. 



r2/4 



A Vr74-hA2 



(11) 



This means that the vector of expectation values 
evolves as (y{t + t)) = hrSriyit)) with L,- = 

diag(Vl-?yr, VI -^r, Vl-e, VI -e)- 

The fraction tJt of atoms that have decayed represents 
a loss of collective squeezing because its correlation with 
the other atoms is lost, whereas it still provides a con- 
tribution ti^/A per atom to the collective spin variance. 
We may use the symmetry of the collective spin operator 
under the exchange of particles to express the mean value 
of, e.g., J2 as (Jf) = f 7V,t + ^N,,{N,, - l){aPai''>) 
where we have used that ((crl )^) = ((ci )^) = 1, and 



(a«ai^)) 



(ai^Vi^)) 



for all i and j (i ^ j). We may 
solve the equation for the correlations between the dif- 
ferent spins 



(aWal^)) 



(J 



—N 



h^ 



X^at(A^at - 1) 



(12) 



l)(ai^V?) 



During a time interval of duration r, rjrN^t atoms de- 
cay by spontaneous emission. This means that (J^) >— > 

U?) = f iV,t(l - 77,) + f iVat(l - Tlr){Nat{l - Vr) ~ 

+ -^rjrN^^, where the last term comes from 
the atoms that have decayed. The correlations given by 
Eq. Ijl2|l are inserted, and for large iVat we find 

(Jf)^(4^) = (l-,70^(jf) + ^(i- (1-^^)2) 



(13) 



where the last line follows in the limit of small atomic de- 
cay, ry,- <C 1. To determine the development of the canon- 
ical variables, we also need the behavior of moments of 



the type {Jx): (Jx) ^ {J'x) = (1 - 'nT){Jx)- Combining 
this result with Eq. lfT3JI . we find 

(Pat)^ (I'at) = (l-^r)(Pat> + 7J7^ (1^) 

and a similar expression for x^t- 

The photons that are absorbed do not contribute to 
the collective Stokes vector, and we find by an analysis 
similar to the above, that 



{S!)^{S'^') = il~ens!) + ih'N,j4)eil-e) 



(l-e)2(5f) + (ft27Vpj4). 



(15) 



in the limit of small e. For the effective Pph variable, we 
find 



{pL)^{pl) = ii-^){plu) + 



(S'x) ' 



(16) 



and a similar expression for Xph- 

Using Eq. (I14II and Eq. Ijl6|l and similar expressions 
for the other elements of the covariance matrix, Eq. Q 
generalizes to 



-f{t + T)=LrSrj{t)S^Lr 



hN,,, 



-Mr 



hN, 



ph 



■N 



{Jx{t)) 2{Sx{t))- 

(17) 
for r]r,e <^ 1 with M,- = diag(r?r,»7T,0,0), and N = 
diag(0, 0, e, e). The factor hNg,t/ {Jx{t)) initially attains 
the value 2, and increases by the factor (1 — ?7r)~^ in 
each time step r. The factor HNph/{2{Sx{t))) is initially 
unity, and is approximately constant in time since the 
light field is continuously renewed by new segments of 
the light beam interacting with the atoms. An exception 
is the optically thick gas discussed below in Sec. IV Bl 

We note that the present accumulation of noise is based 
on the canonical x and p variables entering the covariance 
matrix, and not on the physical spin and Stokes variables 
for the atoms and the photons, respectively. As discussed 
in more detail elsewhere [l^, this introduces difficulties 
in the limit of large atomic decay probabilities. As long 
as the probability for atomic decay is small during the 
process under concern, however, the present approach is 
highly accurate. This is the regime considered in this 
work. 

In the gaussian approximation, the system is fully char- 
acterized by the vector of expectation values (y) and the 
covariance matrix 7. We probe the system by measuring 
the Faraday rotation of the probe field, i.e., by measur- 
ing the field observable a;ph. Since the photon field is an 
integral part of the quantum system, this measurement 
will change the state of the whole system, and in partic- 
ular the covariance matrix of the atoms. We denote the 
covariance matrix by 



7 






(18) 



where the 2x2 sub-matrix A^ is the covariance matrix 
for the variables 1/1 = {xat,Pat)'^, ^j is the 2 x 2 co- 
variance matrix for y2 — {xph,Pph)^ , and C-y is the 2x2 



correlation matrix for yi and y^. An instantaneous mea- 
surement of Xph then transforms A^ as [l^ [la, llJl 



A' 



C-y(7rB^ 



(19) 



where tt = diag(l,0), and where ()~ denotes the Moore- 
Penrose pseudoinverse. 

After the measurement, the field part has disappeared, 
and a new beam segment is incident on the atoms. This 
part of the beam is not yet correlated with the atoms, and 
it is in the oscillator ground state, hence the covariance 
matrix 7 is updated with A^, C^ a 2 x 2 matrix of zeros, 
and BI!y — diag(l,l) before the next application of the 
transformation of Eq. (17). 

Unlike the covariance matrix update, which is indepen- 
dent of the value actually measured in the optical detec- 
tion, the vector (y) of expectation values will change in 
a stochastic manner depending on the outcome of these 
measurements. The outcome of the measurement on a;ph 
after the interaction with the atoms is random, and the 
actual measurement changes the expectation value of all 
other observables due to the correlations represented by 
the covariance matrix. Let x denote the difference be- 
tween the measurement outcome and the expectation 
value of a;ph, i.e., a gaussian random variable with mean 
value zero and variance 1/2. The change of (yi) due to 
the measurement is now given by: 



{yi) 



(yi) = (yi>+Cy(7rB7r)-(x, 



(20) 



where we use that (ttBtt)^ = diag(i3(l, 1)^^,0), and 
hence the second entrance in the vector (x, •) need not 
be specified. 

The gaussian state of the system is propagated in time 
by repeated use of Eq. (17) and the measurement up- 
date formulae 1(1911 - 1(2(1}) . This evolution is readily im- 
plemented numerically, and the expectation value and 
our uncertainty about, e.g., the value of the squeezed p^t 
variable of the atoms are given by the second entrance 
in the vector of expectation values (2/2) = (Pat) and the 
covariance matrix element Aj{2,2) — 2Var(pat)- 

We conclude this section by noting that if one asso- 
ciates with the precise measurement of Xph an infinite 
variance of pat and a total loss of correlations between 
Pat and the other variables due to Heisenberg's uncer- 
tainty relation, the Moore-Penrose pseudoinverse can be 
written as a normal inverse of the covariance matrix, 
(ttBtt)- =diag(B(l,l),oo)-i. Equations (O and ^ 
are then equivalent with the results for the estimation 
of classical gaussian random variables derived, e.g., in 
Ref. fil. 



IV. HOMOGENOUS LIGHT-ATOM COUPLING 

The time evolution of the atomic Pat variable is com- 
pletely determined by the update formulae for the covari- 
ance matrix Ijl7|l and the measurement update formula 



II19II . In the limit of infinitesimal time steps, these for- 
mulae translate into differential equations, and we obtain 
the following equations for the variance of Pat(oc Jz) 

^Var(p,t) = -2k^ (Var(pat))' , (21) 

and 

^Var(pat) = -2K2(l_e)e-''*Var(pat)'-^Var(pat)+?7e''*, 

(22) 
corresponding to the cases where atomic decay and pho- 
ton absorption are neglected and included, respectively. 
Here the light-atom coupling k is given by 



^'-N.,HB' 



(23) 



with X = g^T = ^^^%2 ■ Equation |(2U is readily solved 
by separating the variables, and we obtain 

1 



Var(pat) = 



2K2i + l/Var(pat,o)^ 



(24) 



where Var(pat,o) = 1/2 is the variance of the initial min- 
imum uncertainty state. 



To solve ll22|l . we introduce the change of variable 
Var(pat) = e~''*Var(pat)i and obtain 



— Var(pat) = -2k^(1 - e)Var(pat) - 2?7Var(pat) + V, 



(25) 



which is separable. With (3 = W ^2(i„g) f ^2(i_g) + 2j, 
the solution of Eq. (l22|l reads 



Var(pat) 



P / Var(pat.o) + 5;3f^ + f + e-2^''^(i-)* (Var(pat.o) + j,.^ - f ) ' 



Var(pat,o) 



V 



2K2(i_e) 



/3 _ „-2/3K2(l_e)t 



[Var(pat,o) 



2K2(i_e) 



.vt 



.vt 



2k2(i _ £) 



r 



(26) 



FigureQ]shows the spin squeezing as a function of prob- 
ing time. When atomic decay is not included, the uncer- 
tainty in Pat is a monotonically decreasing function with 
time. When decays are included, a minimum is reached 
whereafter the degree of squeezing starts to decrease. On 
the time scale of the figure, which is chosen to refiect real- 
istic experimental time scales, the increase in Var(pat) is 
hardly visible. From Eq. Ij22|l . we find that the minimum 
in the variance occurs at the instant of time 



in 



^Var(pat,o)+ 2k^~ 



(27) 



In 



4/3k2(i 



Var(p, 



'at.O^ 



2rt^(l-e) 



+ 



/3 



V 



In the typical experimental situation, 77/2^2(1 — e) <^ 
1 which means that (3 ~ \J2t]I{\ — e)/n. In this case 
Eq. lf27jl simplifies to 



1 



2^2^ 



■In 



4^2(1 -£)At 



(28) 



From Eq. l(28|l . we see that t,„in decreases for increasing 
coupHng strength k, and for increasing decay rate t]. In- 
terestingly, the instant of time for the minimum in the 



variance is independent of the initial uncertainty in the 
atomic variable Pat • 

We may now go back to Eq. lf22ll and evaluate the value 
of the variance at time imin- In the regime considered 
above, and in the figure, we find 



Ap(in 




'7 



K V 2(l-e)" 



(29) 



This clearly shows that the higher coupling and the lower 
decay, the better spin squeezing. It is the term linear in 77 
in Eq. lf22|l that is responsible for the 'saturation effect' 
in the variance at early times where the exponential is 
still close to unity, e*** ~ 1. 

To specify for a given number of atoms, how many 
photons we need to obtain optimal spin squeezing in time 
tmin limited perhaps by other experimental constraints, 

we express r\ — *&5s'' *^ ~ '^'i'S^^at, and insert in 
Eq. 1I28I1 . The slow logarithmic dependence and factors 
of order unity can be neglected, and we can introduce e 

1 



via the relation e 



^atfs^, and find $i„i. 



4iVa,. 



If we accept photon absorption at the percent level, we 
obtain 



$innin > 100 




(30) 
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Figure 1: Uncertainty of pat as function of time. The effective coupling is k,^ = 1.83 x 10® s~^. The lower curve is without 
inclusion of atomic decay, and the upper curve includes atomic decay with a rate rj = 1.7577 s~^ and photon absorption with 
e = 0.028. These values correspond for example to a 2 mm interaction area, 2 x 10^ atoms, 5 x 10^ photons s~^, 10 GHz 
detuning, and 852 nm light, appropriate for the ^^^Cs(65'i/2(-F = 4) — 6Pi/2{F — 5)) transition with decay rate 3.1 x 10^ s~^ 
and corresponding atomic dipole moment d = 2.61 x 10~^® Cm. Factors of order unity related to the coupling matrix elements 
among different states of the actual Zeeman substructure are omitted. 



In our case, we have A/ a ~ 1.7 x 10^. A realistic upper 
limit for i,„in is 1 ms, and from Eq. Ij30|l it then follows 
that the photon flux should fulfill 



$ > 10** 



(31) 



then reads 

y = (s^at.ljPat,!: ■ ■ ■ J 2:^at,niPat,n, 2;ph,Pph) ■ (32) 

The generalization of Eq. J^l, to the case with inhomo- 
geneous coupling reads 



V. INHOMOGENEOUS LIGHT-ATOM 
COUPLING 

We now consider two scenarios leading to inhomoge- 
neous light-atom coupling, a case recently discussed the- 
oretically in the literature jlj. First, we shall study the 
case where the coupHng is inhomogeneous as a conse- 
quence of a variation in the intensity of the light beam 
across the sample. Second, we shall consider the case 
of an optically thick sample where the photon field, and 
therefore the coupHng, changes through the atomic sam- 
ple due to absorption. Both cases are readily handled 
within the gaussian approximation. 



A. Case (a): optically thin sample 

We consider the case where the atomic gas is divided 
into, say n, sHces each with local Hght-atom coupHng 
strength Ki . The 2n + 2 column vector of gaussian vari- 
ables describing the 2n collective canonical position and 
momentum variables for the atoms, and the 2 collective 
position and momentum variables for the photon field 



Ht = ^ / f^r,! 



Pat A Pph, 



(33) 



where the summation index covers the different groups 
of atoms. 

To model the effect of an inhomogeneous coupling of 
the light to the atomic sample, we consider n — 10 
different values of k^ chosen uniformly in the interval 
[kI{1-6);kI{1+S)] with S = {0, 0.1, 0.5}. In this way, the 

effective coupling constant \Yl!j=i '*? remains constant 
while the variance in the coupling constants increases. 
The values of the coupHng strength could, e.g., differ be- 
cause of the transverse intensity profile of the laser beam. 
As a consequence, the values of the atomic decay rate 77 
(also proportional to intensity) are different in each sHce. 
The measurement is described by the method in Sec. lIIII 
and the propagation is given by a modification of Eq. Ijl7|l 



7(f + t) = l^s^7^sJ:l^ + m^ + n, 



(34) 



where the (2n -f 2) x (2n + 2) matrix S,- is obtained from 
the time evolution of the system as in Sec. lIIIL and where 



Lr = diag(^l -rir,l, ^1 -VtA, • • ■ , ^/l - ?7r,n, y/l - rir,n , VT 



e,VT~e), 



M.r = hx diag(-Tj — r 
and N = diag(0,0, . 



A^at,!?)! -'Vat,l'7l 



N^t.^V,^ N^t.^^Vr, 



,0,0), 



(J^.l) ' ■ • • ' (J^,„> ' (J. 

.,0, 0,e,e). For convenience, we 
assume that the number of atoms Nat,i subject to a given 
coupHng strength k^ is simply N^^^/n. 

The atomic covariance matrix now has dimension {2n x 
2ri), and it contains the variances of the atomic observ- 
ables in each sUce and the correlations between them. 
Collective observables are described by linear combina- 
tions of the {xat,i,Pat,i) and their variances can be ob- 
tained explicitly. 

From the Hamiltonian ll33ll . it is clear, that the 
probe field couples to the asymmetric collective variable 
12'i=i'^T,iPat,i- The corresponding asymmetric collective 
harmonic oscillator variables involved in the spin squeez- 
ing are, accordingly 



(-'^eff, Pe«) — 



i—l ^i^at^i 2-^i—l ^iPat,i 



VE'Li i^'i ' V^i 



(35) 



The symmetric collective variables that are usually con- 
sidered (see, e.g., the discussion in Ref. [ij and references 
therein), are, on the other hand, given by 



{X,P) 



i=l 



1 " \ 



(36) 



and it is interesting to see how these two sets of variables 
are connected. A straightforward calculation shows that 
we may express the latter variables as 



{X, P) =. a{X,ff, Peff) + b{X^,P^), 



(37) 



where {X± , P± ) are canonical variables which commute 
with (Xeff, Pes) and with the interaction Hamiltonian of 
Eq. (33), and where the coefficients are given by 



V" K 



and 



b{X^,P^) = 




(38) 



, Pat,j)- 



(39) 

From Eq. ll37ll , it follows that the variances of X and P 
may be expressed as 



and 



Var(X) = a2Var(Xeff) + (1 - a^)/2, (40) 



Var(F) = a2Var(PefF) + (1 - a')/2, (41) 



where we have used that 1 = a^ -(- 6^ and that the com- 
ponents {X±,P±) are unaffected by measurements, so 




Figure 2: Uncertainty of the maximally squeezed component 
of the atomic gas (-Pat) as function of time. The higher-lying 
curves show the uncertainty in the symmetric collective pa- 
rameter (Eq. (O) for n — 10 uniformly distributed values of 
K^ in [0.9«:o;l-l«o] (middle) and [0.5«o; I.Skq] (upper) . The 
central effective coupling is kq = 1.83 x 10 s~ , and all other 
parameters are as in Fig. 1. The lower curve is the smallest 
eigenvalue of the covariance matrix, which is the same for the 
two ranges of k^ to the precision visible in the figure. 



Var(Xj^) = Var(PL) = 1/2 for all times (if atomic decay 
is not taken into account). 

In Fig. [21 the lowest curve shows the smallest eigen- 
value of the covariance matrix as a function of time. 
The associated eigenvector represents a combination of 
the canonical variables for the different slices, which is 
maximally squeezed. For the present values of the noise 
parameters (77 and e), we have an overlap very close to 
unity between the eigenvector of this curve and the ef- 
fective asymmetric collective variable P^g of Eq. II35|I . 
This means that this component is indeed the one that 
is maximafiy squeezed. The analytical result for the 
squeezing of this component is obtained from IJ26II with 

For the values for atomic decay and 



photon absorption considered in the figure, also the for- 
mula l|lTll reproduces the fully numerical calculations for 
the symmetric collective coordinate P of Eq. l(3fijl . 



B. Case (b): optically thick gas 

We now turn to the situation where the sample is op- 
tically thick. The probability e for absorption of photons 
through the gas is then larger than, say, a few percent. 
This means that the condition e ^ 1 which was assumed 
in the derivation of the effective light-atom coupling of 



Eq. lOl is no longer fulfilled. By slicing the gas into 
pieces labeled by i = l,2,...,n, within each of which 
the constraint on atomic decay and photon absorption 
rji^ei <ti 1 is fulfilled, we may, however, still locally for 
a fixed slice i use the effective Hamiltonian and address 
the problem in the gaussian approximation. The vector 
of variables describing the system is then of the same 
form as in Eq. (EH, and the Hamiltonian is given by 
Eq. (|SS|l- The considerable absorption of photons from a 
beam segment on its way through the atomic gas, means 
that the update formula for the covariance matrix needs 
to be iterated according to the different local noise and 
coupling strengths. Accordingly, as each beam segment 
passes through the atomic gas for i = 1 to n, we go 
through the following update formulae for the covariance 
matrix (fTzIl : 



7i — Lr,iSr,i7i-lS^ jLr,i 



(42) 



hN, 



(J^At)) 



Mr, 






■N,; 



where the transformation matrix Sr.i is given by a ma- 
trix with off-diagonal elements Kr,i at entrances ((2i — 
l),{2n + 2)) and ((2n-|- 1), 2i)). For example, 8^,2 for the 
case of only two slices (n = 2) is given by 







/ 1 
10 


\ 







S.,2 = 


10 
10 

Kr.2 1 

Vo 


Kt.2 




1 / 


(43) 


The 

must 


constraints on decay 
be fulfilled riT,i,£i '^ 

■■■,^/'^-Vr..i,^/'^-Vr.iA,- 

= diag(0, ...,0,7yr,j,'7 
diag(0,...,0,e„e,)- The 


an( 
1 

r,^,0,. 

full 


1 absorption 

and Lr,i — 


diag(l 

Mt., 


l-e„Vl-e^), 
..,0))., and 
:ovariance ma- 



trix is updated every time the pulse segment passes a 
new slice. When the pulse segment has finally left the 
gas, it is being measured, and 7„ is modified (7„ -^ 7'^) 
according to Eqs. (18) and (19) of Sec. Ill, with the 
2n X 2n submatrix A^ the covariance matrix for the 
variables yi = (a;at,i>Pat,i' • • • '2;at,n,Pat,«)'^, B^.^ the 
2x2 covariance matrix for 1/2 = (xph,Pph)'^ and C-y the 
2ri X 2 correlation matrix for yi and y|^. When we set 
7o(i + t) — 7^(i), we use Eq. lli^ with i = 1 to n to 
describe the interaction with the next beam segment. In 
reality, the light segment corresponding to any practical 
duration r will be much longer than the entire atomic 
sample, and the interaction with one group of atoms has 
not finished before the interaction with the subsequent 
group starts. It is not difficult to see, however, that if 
the atomic dynamics is entirely due to the interaction 
with the optical field, there is no difference between 
the achievements of the real system and those where 
we imagine the atomic sHces separated by free space 
separation distances larger than ct, described precisely 
by the above formulation. In Eq. II42|I . 



For convenience, we give the time and space (slice) de- 
pendence of the parameters in Eqs. lli2|l and l(i3ll explic- 
itly. The change in the classical Stokes vector through 
the different slices due to photon absorption is given by 



\^x.\ 



{Sx,t=o) exp{- ^ ei>), 



(44) 



where the absorption probability in slice i is e^ , and hence 
the total photon absorption probability in the gas is 
(1 — exp(— X]r=i ^«))- The change in {Jx,i) due to atomic 
decay is given by 



(JxAt)) = {JxAo))eM-mt)- 



(45) 



The atomic decay rate rji is a decreasing function of the 
slice-number since fewer and fewer photons are available 
to excite the atoms 



Vt =Vo 



exp(-^Q')> 



(46) 



and finally, the light-atom coupling constant k will de- 
pend on both time and space 



K {t, i) — Koexp(— 



lexp(-?7ji). 



(47) 



j'=i 



where Kg is given as in Eq. ll23ll and every slice contains 
-^at,i = ^at atoms. From the above relations and the ini- 
tial conditions {SxA = ^-^ph,i/2 and {Jx,i) — fiNg_t.i/2 it 
follows that the pre-factors on the noise terms in Eq. Ij42|l 
are given by 



hN, 



\'J X^'i 



-Vi^ ' 



and 



hK 



ph. I 



2{Sx 



=^5ZJ/=i I 



(48) 



(49) 



We have modeled the effect of photon-absorption in- 
duced inhomogeneous light-atom coupling using the pa- 
rameters detailed in the caption of Fig. El The photon ab- 
sorption is varied by varying the detuning, and the Hght- 
atom coupling strength n^ and the atomic decay proba- 
bility 77 are kept constant at the values used in Figs. Q] 
andElby adjusting the photon flux inversely proportional 
to changes in the detuning squared. Figure shows the 
uncertainty of the maximally squeezed component of the 
sample as determined by the smallest eigenvalue of the 
covariance matrix. We see as expected that the degree 
of squeezing decreases with increasing photon absorption 
probability. 

In Fig. 01 we compare for two representative cases from 
Fig- El the uncertainty of the maximally squeezed compo- 
nent of the gas with the uncertainty of the collective inho- 
mogeneous variable P^g of Eq. (J35II . The variance of the 







to the fully numerical result. Only for high photon ab- 
sorption the effects of noise and differences in coupHng 
strength lead to a significant deviation from the numeri- 
cal result. 



VI. PROBING THE DEGREE OF SQUEEZING 

So far, we have not discussed to which extent the max- 
imally squeezed component of the atomic sample will be 
useful and, e.g., set the limit for the precision obtained 
in a measurement of an interesting physical quantity. To 
investigate this point, we follow the work in Ref. [Ol, and 
consider a situation where (i) the sample is spin squeezed 
for a time periode ti (ii) the spin squeezing is stopped, 
and the sample is subject to a spin rotation, and (iii) the 
system is probed, and the rotation angle is determined. 



Figure 3: Uncertainty of the maximally squeezed component 
of the gas (Pat) as function of probing time for varying degrees 
of photon absorption. The percentage of photons absorbed is 
indicated at the solid curves. The gas is decomposed in n 
slices each absorbing e^ = 0.028 of the light intensity. From 
the lower to the upper curve the number of such slices at- 
tains the values n — 1,4,8,13,25, and 50. Other physical 
parameters are as specified in Fig. 1. 



10 



f^ 10 



10" 



"v^^. 




-~-~->__;[;3ll|]| 


■ 75.5% 


(3) 
\ 




V 




(1) 





3 
t [ms] 



A. Noiseless case: Analytical results 

We start by an analysis of the simple case correspond- 
ing to a single atomic sample and a single probe field in 
the noise-less limit. From Sec. IIIIL we have at time ii 



Var(pat(ti)) 



1 



2kHi 



(50) 



where we have used that the atoms are initially in a co- 
herent state with variance 1/2. Since Var(a;at)Var(pat) = 
1/4 in this noiseless case, we also have 



Var(a;at(ii)) = K^ii/2 + l/2. 



(51) 



After the time ii, the light-atom coupling is turned off, 
and the system is subject to a rotation around the y axis, 
described by the interaction HT = —9Jy, where 6 = ujT 
is the small angle of rotation resulting from the action of 
the constant rotation frequency uj in time T, and where 
Jy is the y component of the collective spin operator. 
Making the translation to the effective dimensionless po- 
sition operator as in Eq. Q leads to the Hamiltonian 



HT = —hOaXa 



(52) 



Figure 4: As Fig. Elbut only for n = 4 and n = 50. Curves 
(1) and (3) represent the maximally squeezed component of 
the gas. Curves (2) and (4) display the uncertainty in the in- 
homogeneous collective variable Peff of Eq. H35|l as a function 
of time. 



latter variable can be calculated straightforwardly from 
our knowledge of the time-dependent light-atom coupHng 
constants Hi and the full covariance matrix: Var(Peff) = 

(E^J '^^{t)Ho{t){{p^PJ) - {Pi) {Pj))) / Y.k ^^^W^' Wc SCC 

that for low and moderate photon absorption, the result 
for the effective asymmetric variable of Eq. Ii;-i5|l is close 



-^. To obtain an estimate for 



where a = y/{Jx)/h - 
the unknown classical variable 9, we follow the ideas in- 
troduced in Ref. Q, and treat the rotation variable 9 as 
a quantum variable within our gaussian description. The 
total system is then described by two atomic variables 
and one rotation variable y = {9,Xat,Pat)'^ ■ The corre- 
sponding transformation matrix follows from Heisenbergs 
equations of motion with the Hamiltonian in Eq. Ij52|l and 
in the basis (6',Xat,Pat) we obtain 



(53) 




10 



from which we verify that, e.g, Pat -^ Pat + o;9. Equation 
(jnj now determines the time-evolution of the system, and 



we find the following covariance matrix at time ^2 after 
the rotation: 



l(t2 



2Var(0o) 









a2Var(0o) 





2Var(2;at 


(ii)) 







a2\ar(9o) 







2Var(pat(ii))+a'2Var(0o 



(54) 



where Var(pat(ii)) and Var(a;at(ii)) are given by 
Eqs. H5()ll and 115 l|l . respectively. 

Finally, at times t > t2, the rotation is turned off, 
and the sample is probed by the light beam as in the 
time interval [0;ii]. The transformation matrix is deter- 
mined by Heisenberg's equations of motion for the vari- 
ables y = {9, Xat, Pat, Xph, Pph)'^ with the Hamiltonian © 
and is given by 



S = 



/ 1 \ 

10 

K 1 

Vo o/ 



(55) 



The covariance matrix of the system is propagated ac- 



cording to Eq. Q. The measurements on the photon field 
are described as in Eqs. Ijl8|l and IjliJII (see also Ref. (Jl). 
The submatrix A^ is now the 3x3 matrix pertaining to 
the variables (6',Xat,Pat), and Cj is the 3x2 covariance 
submatrix describing the coherences and correlations be- 
tween these three variables and the photon field. We are 
interested in the uncertainty on the value of 9, i.e., the 
(1,1) entrance in the covariance matrix. To find this as a 
function of time, we follow the procedure in Sec. lIIII and 



calculate the difference between A^ after n and n-l-1 iter- 
ations, and consider the limit of infinitesimal time steps. 
In general, the differential equations obtained in this way 
are matrix Ricatti equations and may be solved in stan- 
dard ways [2^. In the present case, the solution reads 
for probing times t >t2: 



J 



Var(6l(t)) = Var(6lo) 



Re{{{9-{9)){pat-{Pat)))t.) 
Var(p(t2)) 



1- 



{l + 2Yaiip{t2))K^t~t2)) 



(56) 



r 



where the covariances at time ^2 are given in Eq. II54II . 
We see from Eq. (I5(ill that the variance of the variable 9 
does not decrease forever. In the long time Hmit, we find 



Var(6'(t -> oo)) = Var(0o) 



Var(pat(ii)) 



Var(pat(ii))+«2Var(0o) 



(57) 
This shows, as expected, that the limiting value only de- 
pends on the squeezing and the rotation until time t2. 
For large a parameter (many atoms) and for a sufficiently 
large initial variance of 9, the result in Eq. Ij57|l reduces 
to 



Var(6l) 



Var(pat(ii)) 



(58) 



The ratio of the variances of in a measurement with (S) 
and without (NS) spin squeezing is given by 



S = ^v„(A«.,) 



(59) 



Since Var(pat(^i)) g]0; 1/2] this shows that one may gain 



a significant factor in precision on the variable 9 by pre- 
squeezing the sample. 

Finally, we note that the result of Eq. lf57jl 
may be obtained directly by considering the cor- 
responding classical gaussian probability distribution 
PiPat,0) (X exp (-p2j(2Var(pat)) - 0V(2Var(^))) . As 
a consequence of the rotation, pat transforms accord- 
ing to Pat — > Pat + a^j and therefore the prob- 
ability distribution after rotation reads P(pat,^) c< 
exp {-{pat - a0)V(2Var(pat)) - 0V(2Var(0))). A mea- 
surement of the variable Pat leads to a distribution in 9 
only, from which the variance of 9 is read off with the 
result given in Eq. H57II . 



B. Noise included: Numerical results 



Whereas in Sec. IVI Al it is clear that it is the col- 
lective variable Pat that is squeezed, in the case of an 
atomic ensemble with an inhomogeneous light-atom cou- 
pling we only know from the analysis of Sees. IV Al and 
IV Bl that there exists a component that is squeezed, and 
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that this component for moderate noise is very accurately 
approximated by the asymmetric collective variable Peff 
of Eq. ll35ll . The question we address now is whether 
it is the variance of this component that will show up 
in a measurement of a classical parameter, such as the 
rotation parameter 6. 

The formalism necessary for handling this problem was 
developed in Secs. lIIII and IV Bl In short, for n slices of gas 
each fulfilling 6^,77^ <C 1, (« = 1, . . . , n) we first propagate 
and perform measurements on the system of 2n collective 
atomic position and momentum variables and 2 collective 
photon position and momentum variables. At time ti, 
the light field is turned off, and the atomic sample is 
for t e [ti]t2\ subject to a rotation around the y axis 
described by the effective Hamiltonian 



Ht =^ -ne^aiXs^t.i, 



(60) 



with 6 — ijjT as in Sec. IVI Al and with coupling constants 
ai determined by a generalization of the result in Eq. lf52|l 



\^x^\ 



N,, 



-riiti 



(61) 



Spontaneous emission of photons is neglected in our ap- 
proach, so the Qfi's are fixed by their values at the 
instant of time ti when the photon field is switched 
off, and the possibility for stimulated atomic decay 
disappears. The transformation matrix S correspond- 
ing to the Hamiltonian in Eq. Ij6()|l is readily found 
from Heisenberg's equations of motion for the variables 
(^,a;at,i'-Pat,i>-- •'2;at,n,Pat,n)- Its diagonal entries are 
unity, the i = 1, . . . , n {0,Pat,i) entries are assigned the 
values Ui, and the rest are zero; a natural generaliza- 
tion of Eq. lf53ll . The propagation in time of the co- 
variance matrix is then determined by Eq. Jsj. At time 
t2 the rotation is stopped, and for times t > t2, the 
atom-light Hamiltonian is turned on again. First the 
initial covariance matrix for 6, atomic slices and the pho- 
ton field (0,Xat,l,Pat,l,-- •,a;at,n,Pat,n,a;ph,Pph) IS Set Up. 

This involves the covariance from the previous part sup- 
plemented by the position and momentum variables of 
the photon field. The dynamics of this enlarged covari- 
ance matrix is described by suitable modified versions of 
Eqs. (17) and (19) of Sec. HI. 

propagation with this enlarged covariance matrix is 
performed by the standard equation of Sec. lIIII nronerlv 
adjusting the transformation matrix and the matrices as- 
sociated with noise. 

We aim to extract from our numerical study that the 
variable of relevance in the probing of the rotation an- 
gle is the maximally squeezed component, i.e., at mod- 
erate noise levels it is essentially the optimally squeezed 
asymmetric Peff variable of Eq. l(35|l and not the symmet- 
ric collective variable P of Eq. l(36)l . From Heisenberg's 
equation of motion it follows that P^g and P transform 
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Figure 5: Uncertainty of the parameter S as a function of time 
and for different variances of the coupling strength as specified 
in the text. The gas is sliced in n = 10 pieces. The number 
of atoms and photons are as in the preceding figures. The 
value of ai is 0.2236. The dashed curves show the limiting 
uncertainty in the 9 parameter as estimated from Eq. I|65|l for 
the standard collective variable P of Eq. IHfill with the small- 
est variance in the coupling strength for the lowest dashed 
curve and the highest variance for the upper dashed curve. 
The constant horizontal solid line gives the limiting value of 
Eq. (64), as obtained by the maximally squeezed component 
Peff of Eq. IjnSJ, and it is independent of the variance of the 
coupling strength. The decreasing full curve is a collection of 
indistinguishable curves showing the numerical results for all 
the different variances of the coupling strength (see text). 



according to 



eff -^ -< eff 




(62) 



and 



P ^ P 



En ' 



(63) 



A generalization of the result in Eq. H58II then yields the 
following expressions for the variance of in the long time 
limit 



Var(6l) = 



and 



Var(6l') 



Var(Peff) 



Vsi{P) 



(64) 



where Var(P) is given by Eq. Ij41|l . 



(65) 
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Figure [3 shows results for inhomogeneous coupling 
modeled by choosing n = 10 different values of k^ uni- 
formly over the interval [(1 — (5)K§/n; (1 + (5)Kg/n] with 
6 G {0, 0.02, 0.1, 0.2, 0.3, 0.4, 0.5}. As in Sec.El the effec- 
tive coupling strength is fixed by kq. In the figure, the 
soHd lines are independent of fluctuations in the coupHng 
strength. The lowest solid fine shows the asymptotic un- 
certainty of 6 as obtained by Eq. Ij64|l . The decreasing 
solid curve is the numerical result, converging towards 
this value. It represents a collection of indistinguishable 
curves showing the numerical results for all the different 
variances of the coupling strength. We observe that the 
decreasing solid curves show a better estimation of the 
rotation angle 9 than the prediction by the symmetric 
collective variable shown by the dashed curves in the flg- 
ure. The fact that the decreasing full curves converge to 
the value determined by the maximally squeezed com- 
ponent signifles that this indeed sets the limit for the 
precision of the measurement. 



VII. CONCLUSIONS 



The gaussian approximation for the collective quan- 
tum parameters including possibly an external classical 
parameter allows us to include the measurement process 
directly, and to obtain analytical results in the noise-less 
case and in the limit of low noise. Also the theory is read- 
ily generaHzed to handle situations which have resisted 
a satisfactory treatment with other theoretical methods. 
For example, the case of an optically thick gas with cor- 
responding inhomogeneous light-atom coupling can be 
treated and even understood analytically to a large ex- 
tent. 

We have shown that in the present case of squeezing 
of the spin of an atomic ensemble by using a continu- 
ous wave coherent light beam, it is indeed the maximally 
squeezed component of the atomic gas that determines 
the precision with which one can estimate the value of 
an external perturbation. 

At present, we seek to address a series of other prob- 
lems in continuous variable quantum physics including 
generation and detection of flnite band-width squeezed 
light and estimation of time-varying external perturba- 
tions. 



In this work we have given a comprehensive account 
of the theory of probing and measurements in the gaus- 
sian state approximation. We have followed the ideas of 
Refs. (2,llJ|i ^'^d we have provided a complete analysis of 
the method and its strengths by analyzing in detail the 
problem of spin squeezing. 
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